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Abstract 

When subject to a weak magnetic impurity distribution, the order parameter and 
quasi-particle energy gap of a weakly disordered bulk s-wave superconductor are sup- 
pressed. In the Born scattering limit, recent investigations have shown that 'optimal 
fluctuations ' of the random impurity potential can lead to the nucleation of 'domains ' 
of localised states within the gap region predicted by the conventional Abrikosov-Gor'kov 
mean-field theory, rendering the superconducting system gapless at any finite impurity 
concentration. By implementing a field theoretic scheme tailored to the weakly disordered 
system, the aim of the present paper is to extend this analysis to the consideration of 
magnetic impurities in the unitarity scattering limit. This investigation reveals that the 
qualitative behaviour is maintained while the density of states exhibits a rich structure. 

1 Introduction 

In the absence of Coulomb interaction effects, the spectral and transport properties of a 
bulk singlet s-wave superconductor are largely insensitive to the presence of a weak non- 
magnetic impurity potential. This effect, which is ascribed to the Anderson theorem Q], 
limits the influence of long-range phase coherence phenomena to situations in which low- 
energy quasi-particles persist: notably, the physics of hybrid SN-compounds, and those which 
exhibit unconventional (e.g. d-wave) symmetry. However, another method of inducing low- 
energy quasi-particle states in the disordered superconducting environment is to impose an 
external time-reversal symmetry breaking perturbation which has a pair-breaking effect on 
the condensate. 

Considering mechanisms of time-reversal symmetry breaking, it is possible to conceive of 
at least two distinct physical situations. The first is the imposition of a homogeneous magnetic 
field: the diamagnetic properties of the superconductor limit considerations to either a vortex 
phase of a type II superconductor, or to superconductors whose lateral dimension is smaller 
than the penetration depth. In each case, field lines are able to penetrate the sample. A 
second method of breaking the intrinsic time-reversal symmetry of the system is to impose a 
magnetic impurity distribution. With both mechanisms, the perturbation acts with opposite 
sign on the two members of the Cooper pair. In the first case, the paramagnetic term in 
the single particle Hamiltonian reverses sign under p — > — p. In the second case, the spin 
of the magnetic impurity acts on different spin components with different sign. Building on 
the existing literature, the aim of this paper is to explore the quasi-particle properties of a 
weakly disordered superconductor subject to a magnetic impurity distribution with unitarity 
limit scattering. 

1.1 Background: Abrikosov-Gor'kov Theory 

In the earliest work in this area, attention was focussed on the infiuence of a weak magnetic 
impurity distribution in which the influence of disorder could be treated in the Born approxi- 



mation. In a seminal work by Abrikosov and Gor'kov it was shown that the pair-breaking 
potential brings about only a gradual suppression of the superconducting order parameter. 
More surprisingly, according to the self-consistent mean-field theory, the energy gap in the 
quasi-particle density of states (DoS) is suppressed more rapidly than the order parameter, 
admitting a region in the phase diagram where the superconductor exhibits a 'gapless' phase. 
More precisely, defining the dimensionless control parameter, 

^ r,|A| ' 

where |A| represents the self-consistent bulk order parameter, and I/ts denotes the Born 
scattering rate due to the magnetic impurities, the Abrikosov-Gor'kov mean-field theory 
shows the energy gap to vary as eg = |A|(1 — (^2/3^3/2^ showing an onset of the gapless region 
when C = 1- Soon after its introduction, it was realised that the general Abrikosov-Gor'kov 
scheme applies equally to other mechanisms of pair-breaking (such as that imposed by a 
uniform magnetic field in a thin film or by a supercurrent) — requiring only a reinterpretation 
of the dimensionless parameter C (see, e.g., Ref. [||). 

In later works ||, ||, ^, ||] , various authors explored the influence of isolated magnetic 
impurities. In particular, in the unitarity limit, it was shown that a single classical magnetic 
impurity of spin S leads to the local suppression of the order parameter |Q, || and nucleates 
a bound sub-gap quasi-particle state at energy 

Cb _ |1 - a| 

where, defining the density of states 1/{L'^S) of a normal conductor, with S being 

the single-particle level spacing, a — (7ri/JL'*|S|)^ represents the dimensionless scattering 
amplitude associated with the magnetic impurity. 

For a finite impurity concentration, the sub-gap states weakly overlap, hybridise, and 
broaden into a band ||] centered on energy Cb- Here, as in the Abrikosov-Gor'kov theory, the 
mean- field theory again predicts a gradual suppression of the quasi-classical energy gap e^, 
with the superconductor entering the gapless phase when = (^q = (1 — a)^ . Increasing the 
magnetic impurity concentration a, two qualitatively different situations can be realised: in 
the first case, when C reaches the value ^1 , the impurity band can merge with the continuum 
of bulk states before the system enters the gapless phase (Ci < Co), while, in the second case 
(Co < Ci), the opposite situation pertains (see Fig. [l]). 

1.2 Beyond Mean-Field Theory 

Despite the success of the Abrikosov-Gor'kov mean-field theory and its extension to the 
unitarity limit scattering, two questions present themselves: 

• Firstly, the existence of phase coherent low-energy quasi-particle states in the gapless 
phase renders the spectral and (thermal) transport properties of the superconductor 
susceptible to the influence of long-range quantum interference effects. 

• Secondly, according to the mean-field description, a hard energy gap is maintained 
up to a critical concentration of magnetic impuritiesn. Yet, being unprotected by the 



^ In the Born scattering limit at T = 0, the critical concentration of magnetic impurities at which the 
energy gap goes to zero is 2e~'^/'* ~ 0.91 times the critical concentration at which superconductivity is 

destroyed. In the unitarity limit, the critical value is decreased to 2[(1 — a)/(l + a)] e 1+° . 
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Figure 1: Quasi-particle density of states of a disordered superconductor with magnetic impurities 
drawn from a Poissonian distribution; (a) a — 0.25 (^i < ^o), the value of is increased from = 0.05 
to ( — 0.6. The localised excited state for the one impurity problem is located at Eb/IAI — 0.6. (b) 
a = 0.6 ((^0 < Ci)5 the value of is increased from — 0.05 to = 0.9. The localised excited state 
for the one impurity problem is located at Eb/IAI = 0.25. 



Anderson theorem, it would seem that the gap structure predicted by the mean-field 
theory is untenable and may be destroyed by 'optimal fluctuations' of the random 
impurity potential. 

In recent years, both issues have come under scrutiny. In particular, it has been shown 
that the long-range, low-energy spectral and transport properties of a weakly disordered 
superconductor can be classified according to the their fundamental symmetries In the 
quasi-classical limit epT 3> 1, where 1/r represents the scattering rate of the non- magnetic 
impurity distribution, the spectral and transport properties of the weakly disordered system 
can be presented in the framework of a statistical field theory of non-linear a-model type (for 
a review see, e.g. [0). Within this approach, the conventional Abrikosov-Gor'kov mean- 
field theory is identified as the set of (homogeneous) saddle-point equations. Mesoscopic 
fluctuations due to quantum interference effects in the particle/hole channel are recorded in 
the soft field fluctuations around the homogeneous mean-field solution. 

In the present case, such investigations reveal that mechanisms of quantum interference 
lead to the delocalisation of the quasi-particle states even in low dimension |l^ |l^ . 
This behaviour provides a striking contrast with that of other superconducting (and nor- 
mal metallic) systems where mechanism of quantum interference have a tendency to bring 
about localisation of the quasi-particle states. At the same time, the same general theoretical 
framework provides a means to explore the integrity of the gapped phase in the supercon- 
ducting system, 'optimal fiuctuations' of the random impurity potential(s) nucleate domains 
or droplets of localised tail states below the predicted mean-field gap edge Such tail 

states are accommodated by instanton field configurations of the non-linear cr-model. 

The tail states predicted by the quantum field theory differ substantially in character 
from the bound states induced by an isolated magnetic impurity. The former derive from 
mesoscopic fluctuations of the electron and hole wavefunctions of the normal system: specif- 
ically, in regions where the phase sensitivity is anomalously high, the pair-breaking effect 
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of magnetic impurities (or an external magnetic field) is enhanced over that predicted by 
the mean-field. Here, quasi-particle states localise over length scales comparable to the su- 
perconducting coherence length, ^ = (D/|A|)^/^, where D = v^r/d represents the classical 
diffusion constant. 

As well as presenting a concise review of the field theory of the disordered superconducting 
system with magnetic impurities, the aim of the present paper is to explore the integrity of 
the sub-gap state picture when in the presence of unitarity limit scattering. More precisely, 
at the level of mean- field, we have seen that, over a wide region of the phase diagram, the 
latter induces a delocalised band of bulk states centered on the bound state energy e^. In 
this case, do optimal fluctuations of the random potential lead to the nucleation of localised 
states in the vicinity of the narrow band? 

The paper is organised as follows. In section || we refine the field theory of the weakly 
disordered superconducting system to incorporate the presence of magnetic impurities. Here 
we will take the magnetic impurities to be drawn from a random Poisson distribution al- 
lowing a continuous interpolation from the unitarity scattering limit to the Born scattering 
limit. Having obtained the low-energy effective field theory, in section 2.6 we will explore 
the homogeneous mean-field solution obtained from the saddle-point of the effective action. 
In this case, we correctly recover the phenomenology of Ref. and identify the limit in 
which the Born scattering Abrikosov-Gor'kov theory |^ is obtained. By exploring instan- 
ton field configurations of the action, in section || we explore the integrity of the mean-field 
density of states. In particular, we will show that the gap edges predicted by the mean-field 
theory become mobility edges separating regions of bulk delocalised states from localised 
'droplet' states generated by optimal configurations of the random impurity potential. A 
brief discussion of these results is contained within the concluding section. 



2 Field Theory of the Superconducting System 

Previous investigations have shown that, in the quasi-classical limit, the properties of the 
weakly disordered superconducting system can be expressed in the framework of a statistical 
field theory of non-linear cr-model type. In the present case, one must consider simply 
how to tailor this analysis to the consideration of unitarity scattering limit of the magnetic 
impurity system. Since the general theoretical framework has been reviewed in a number 
of publications |ll| and discussed for the magnetic impurity system in particular p5| ], 
we will keep our discussion concise focussing primarily on the idiosyncrasies of the present 
theory. 



2.1 The Model 

In the mean-field BCS approximation, a bulk s-wave disordered superconductor in the pres- 
ence of magnetic impurities is specified by the Gor'kov Hamiltonian 



H |A|ar 



where the index PH refers to the particle/hole space, and Pauli matrices a^^ operate in the 
spin space. Here 

^ = Cp + ^(r) + ^S(r) • cr--^" . 

denotes the single particle Hamiltonian, with = p^/2m — ej?, ep is the Fermi energy, 
and |A| is the spatially homogeneous order parameter determined from the self-consistency 
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condition, A = —L'^gA{ipiipi)- In addition to a non-magnetic impurity potential, V{v) drawn 
at random from a Gaussian white-noise impurity distribution with zero mean and variance, 

{V{r,)V{r2})v = -^S{r,-r2) , 

the system is subjected to a classical quenched Poisson distributed magnetic impurity po- 
tential: 

S(r) = L^^5'^(r-r,)S, , 

i 

i.e. where the points are drawn from a random Poissonian distribution. Here, for simplicity, 
we will suppose that the spins corresponding to different magnetic impurities are statistically 
independent, and that the distribution over the orientation is uniform, while the magnitude, 
S is fixed: i.e. P({SJ) = Ui^^l ~ S^)- 

Before proceeding, we should comment on the limitations of the present scheme. In 
the unitarity limit, one would expect spin 1/2 quantum magnetic impurities to be fully 
Kondo screened by the itinerant electron system. Our classical model is therefore limited 
to situations in which either the magnitude of the spin is sufficiently large that the moment 
can not be fully compensated or, more realistically, to systems where the mutual RKKY 
interaction of the magnetic impurities lead to a spin glass ordering of the moments. 



2.2 Generating Functional 

To formulate a field theory of the non-interacting superconducting system, we will follow the 
standard scheme ||l^, and begin with the generating functional for the single quasi-particle 
Green function: 

Z[j] = J D{'ip\'i{}) exp |i J dr (^e+ - iJc) + V'^J + J V 

Here e+ = e + iO, while and V'(r) represent two independent eight component super- 

vector fields with ph, SP and boson/fermion (bf) internal indices. By incorporating an equal 
number of fermionic and bosonic components, the normalisation 2[0] — 1 is automatically 
imposed. In the following, it is convenient to implement a gauge transformation ip i-^ Uip, 
where U = E™ - <^ and E™ = diag(l,0)pH, = diag(0, 1)ph, whereupon, the 



Gor'kov Hamiltonian takes the simpler form: 



Cp + V{r) af + lA^r + JS(r) • a' 



As mentioned above, to determine the influence of quantum interference effects on the 
disordered superconducting system, it is useful to first classify the microscopic Hamiltonian 
according to its fundamental symmetries. In the absence of magnetic impurities the Gor'kov 
Hamiltonian exhibits both time-reversal symmetry, and the particle/hole symmetry 



Ha = -a™ . 



^2 



(1) 



In presence of magnetic impurities the time- reversal is broken, while the particle/hole sym- 
metry is conserved. Applied to the corresponding Gor'kov Green function, Gg *(e) — (e± — 
Hg)~^, the particle/hole transformation (||) converts an advanced function into a retarded 
one: 
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As usual ||T^, |ri|, to accommodate quantum interference effects in the particle/hole chan- 
nel, it is convenient to affect a further space doubling tp^G^.{e)'ip — ^G^.{ecr^^'^)'^ , where, in 
the charge conjugation space, the vector fields take the formn: 



V2 ® v^^^ 



This completes the formulation of the generating functional for the single particle properties 
of the Gor'kov Hamiltonian. The theory is specified in terms of 16-component supervector 
fields ^' and ^ with the following symmetry relations. 



with 7 = ia2^Efl - crf'^-BH- As before Efl = diag(l,0)BF and E22 = diag(0, 1)bf represent 
projection operators on the boson/fcrmion space. 

2.3 (Non-Magnetic) Impurity Averaging 

An ensemble average of the generating functional over the non-magnetic impurity distribution 
V induces a quartic interaction of the fields which can be decoupled by means of a Hubbard- 
Stratonovich transformation with the introduction of 16 x 16 supermatrix fields Q(r) 



exp 



DQexp 



dr { ^ str Q2 _ ^^Q^PH^ 



The symmetry properties of Q are inherited from the dyadic product cr™^'(r) (g) 4'(r) and 
impose the condition 



(2) 



In principle we could immediately subject the generating functional to a further average 
over the Poisson distributed magnetic impurity potential. However, such an approach proves 
to be unprofitable. Since the typical separation of magnetic impurities, £s — vpTg, is greatly 
in excess of the mean-free path associated with the non-magnetic impurities, £ = vpT, it is 
more sensible to postpone the second ensemble average until the quasi-classical theory has 
been developed. Therefore, at this stage, let us proceed by integrating out the superfields 
after which the generating functional assumes the form 



J DQexpl^^ J dr strQ^r)-^ J dr str(r|lng-»| 



(3) 



where Q represents the supermatrix Green function, 



JS • cr" 



Go 



^ The transposition operation for the supervectors 1/) and Jp^ and the supermatrix F is chosen according 
to the convention: 
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2.4 Intermediate Energy Saddle-Point: the Non-Linear a-Model 

To make further progress it is necessary to employ a saddle-point approximation. Following 
Ref. ||l6|| , it is convenient to implement a two-step procedure making use of the hierarchy of 
energy scales which place the superconductor in the quasi-classical and dirty limits: 

CF > - > {-,|A|} > S . 

T Ts 

To implement the quasi-classical approximation, we therefore temporarily suspend the en- 
ergy source, e, the order parameter, |A| and the magnetic impurity potential, JS, and seek 
an intermediate energy scale saddle-point. Such an analysis is discussed in detail in the lit- 



erature 1 17 and here we only recapitulate the results. A variation of the action with respect 



to Q obtains the saddle-point equation 

Q(r) = — a3™(r|^o|r). 

Taking into account the analytical properties of the Green function, in the pole approxi- 
mation, one obtains the conventional saddle-point solution Qsp = ® a^'^'. However, the 
saddle-point solution is not unique but spans the non-linear manifold — I: in the absence 
of the external symmetry breaking perturbations, the saddle-point equation admits an entire 
manifold of homogeneous solutions parameterised by transformations Q = TQspT~^, where 
T is a supermatrix, constant in space, and compatible with the charge conjugation symmetry 
properties of Q, (||). 

Fluctuations of Q transverse to the saddle-point manifold are massive and may be inte- 
grated out within the saddle-point approximation, justified by the large parameter 1/t5. By 
contrast, fluctuations which preserve the non-linear constraint = I are massless and must 
be integrated exactly. In the same quasi-classical approximation epT 3> 1, the matrix Green 
function takes the form, 

PH^ /^ri -f-rz 



f/o(i"i,i"2) = -iTTiyfdQri - r2\)(T™Q 

where the Friedel function fd{r) = (9G*(r, 0))y / ('^Gq{0, 0))y denotes the impurity averaged 
single-particle Green function. 

To determine the intermediate energy scale action, we now restore the symmetry breaking 
parameters e, |A| and JS. Expanding the term strln^^^ appearing in (|^), to leading order 
in e and |A|, one obtains: 

Ing-i ~ In^o-i + [^0 (e^f - l^WV")] + In f I - JS • a'A , (4) 



where, importantly, the magnitude of the magnetic impurity potential has been left unre- 
stricted and where the neglected terms turn out to be of order er <C 1 or |A|t <C 1. Taking 
into account slow fluctuations Q{r) — T~^(r)CT™ (g) a^'^T{r), a gradient expansion of the flrst 
two terms of the series (jij) to leading order in |A|t and er recovers the familiar non-linear 
(T-model action for the disordered superconductor, 

Sa[Q] = -Y J str [D{VQf + 4 (ze+a™ ® + |A|a™) Q] , (5) 

where D = Vpr/d represents the classical diffusion constant associated with the non-magnetic 
impurities. Taken together with the spin scattering contribution, the average generating 
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functional assumes the form 
where 

Ss[Q] ^IJdr str(r| In (l - ^o^S • a''') |r> . 

To make sense of the spin scattering component of the action, it is now necessary to implement 
the magnetic impurity average. 

2.5 Magnetic Impurity Averaging 

If we restricted ourselves to the limit of Born scattering, we would be free to expand the 
action S's[Q] to second order^ in S. In this case, once ensemble averaged over the magnetic 
impurity distribution, we would obtain the weak coupling action S[Q] ~ So[Q] + Sg^[Q], 
where 

With the identification (d„ ~ 3) 

TTV , 

2ans = — , (6) 

where l/r^ represents the magnetic impurity scattering rate, a = {ni'JL'^S)^ , and Ug is 
the magnetic impurity concentration, this result coincides with that obtained by Ref. |l5[] . 
However, in the unitarity scattering limit, we are not at liberty to freely take the impurity 
potential from underneath the logarithm. 

To simplify the action Ss [Q] , let us recall that the typical separation of magnetic impuri- 
ties, is, is greatly in excess of the mean- free path associated with the non- magnetic impurities, 
£. (In the opposite limit, superconductivity would, in any case, be fully suppressed.) In this 
case, we may affect the approximation 

gQ{r.,,rj) ~ -i7rw™(3(r.,)% . 

Making use of this relation, an expansion and re-summation of the logarithm leads to the 
result: 

SsiQ] -^Y. '^t^l^ P + (i7rz.JL'^)a™g(r,)S, • a'^] . 

i 

In this approximation, the functional has become separable in the individual magnetic 
impurity scatterers. In principle, we could proceed directly by subjecting the generating 
functional to an ensemble average over the random Poissonian distributed magnetic impurity 
distribution. However, for a general supcrmatrix Q, without expanding the logarithm, the 
average over the independent spin degrees of freedom is laborious and not illuminating. 
Instead, we will follow a different program. 



^ Note that the term Unear in S generated by the expansion vanishes when projected onto the singlet 
saddle-point {pT), considered below. 
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Firstly, by taking into account the symmetry of the soft degrees of freedom, we can further 
simpHfy the analysis: specifically, from the previous analysis of the Born scattering theory, 
it is evident that fluctuations of Q which are not proportional to F'' will be rendered massive 
by the magnetic impurity potential. Since we are interested in the low-energy content of the 
theory, we may therefore specialise our considerations to the singlet degrees of freedom of Q, 
i.e. Q ^ Q ^Y'^ . In this case, if we undertake the trace over the spin degrees of freedom of 
Q, one obtains 

i 

where = {nvJ L'^\Si\Y ^ constant, is the dimensionless scattering amplitude associated 
with each magnetic impurity, and the notation strg indicates that we have carried out the 
trace over the spin indices. 

Secondly, the structure of the non-linear (T-model action (H) shows that the dominant 
contributions to the generating functional arise from field configurations which vary slowly 
at the scale of the coherence length. Therefore, in the dense magnetic impurity limit, where 
the typical separation between magnetic impurities, is — vpTs, is much smaller than ^, one 
can expect the density of spin scatterers, p{r) — (5(r — r^), to be dominated by it's smooth 
average^, Ug. In this approximation, one obtains 

Ss[Q] = Y J rfrstr8ln[I-aafQ(r)afg(r)] . 

Notice that the invariance properties of Q on the saddle-point manifold imply a symmetry 
of the action under the transformation a i— > 1/a: the weak (or Born) scattering limit is 
therefore 'dual' to the strong scattering limit. 

This concludes the construction of the field theory of the disordered superconducting 
system. Single quasi-particle properties of the superconductor are presented as a functional 
field integral (. . . )q = J DQ . . . e~'^'="['31 involving the effective non-linear cr-model action 

Ses[Q]=So[Q] + Ss[Q] . 
In particular, the local quasi-particle DoS can be obtained from the functional integral 

(Ke, r))v,s = ^5R(str8 [^3^ ^ cTg™ ^ arQ{r)])Q . 

Although the soft mode action is stabilised by the large parameter epT ^ 1, the major- 
ity of field fluctuations of the action are rendered massive: the energy source e, the order 
parameter |A| and the magnetic impurity potential all lower the symmetry of the interme- 
diate energy scale theory. To assimilate the effect of these terms, and to establish contact 
with the Abrikosov-Gor'kov theory Q in the Born approximation limit, and with Ref. 1^ 
in the unitarity limit, in the following section, we will proceed by exploring the low-energy 
saddle-point structure of the theory. 



* Such an approximation circumvents the need to implement the averaging over the Poisson distribution 
explicitly. Later, we will see that this brings with it considerable simplification in the instanton analysis, 
while leaving the mean-field analysis unchanged. 
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2.6 Low-Energy Saddle-Point 

To explore the low-energy structure of the theory we proceed by varying the action S'cffiQ] 
with respect to Q subject to the non-linear constraint ~ 1. In doing so, one obtains the 
saddle-point equation 



where the scattering rate 1 /r^ is related to the magnetic impurity concentration Ug and the 
scattering amplitude a through the relation (|^). The latter must be supplemented by the 
self-consistency condition 

— l^l = 5?AE«tr8 [al-®afQ] 
gA opo — 

n 

where is the BCS coupling constant and Q represents the solution of the mean-field 
equation taken at the Matsubara frequencies e t—^ ien = i(2n + 1)-k / [3. Applying the Ansatz 
for the saddle-point, 

Q(r) = sin0"(r)f7™ + cos ^"(r)^!," (g) af' , (7) 

where the matrix 9 = diag(i0BB, ^ff)bf is diagonal in the superspace and independent of the 
other indices, one obtains the saddle-point equation 

T / ' , . , ~\ 1 sin 29 , . 

i:>V26' + 2 lie sin 6*- A cos 61 r = 0. (8) 

V / r,. I + a2 + 2a cos 26* 

As expected, in the limit a <^ 1, an expansion of the spin scattering term obtains the 
Born scattering equation analysed in Ref. Here, taking the saddle-point solution to be 
homogeneous in space and symmetric in the superspace, i^BB — ^ff = 9, one recovers the 
Abrikosov-Gor'kov mean-field equations. Similarly, in the strong scattering limit a ^ 1, the 
action coincides with that obtained in the Born scattering limit. This coincidence reflects the 
duality seen on the level of the action and can be understood qualitatively in the following 
way: when the local spin scattering potential is very strong, the wavefunction for the low- 
energy quasi-particles states is strongly suppressed at the impurity centre. As a result, the 
matrix element for the effect spin scattering rate is itself strongly suppressed. However, our 
main interest is in the crossover region where the dimensionless scattering rate is a ~ 1. Here 
the saddle-point equation is given by (H). 

As with the Born scattering limit, to analyse the saddle-point equations, we first seek a 
homogeneous supersymmetric solution, 0mf = ^mfI^^- Then, defining a 'renormalised' energy 
and order parameter, 

i cos 0MF 

(9) 



lAI = lAI 



2t, 1 

1 sin 



2t, 1 



the homogeneous saddle-point equation ^ takes a conventional BCS form, iesin^MF 
|A| 

cos^MF- Defining the dimensionless parameter u — £/|A| such that 



cosWuiF = , smc/.iF = 
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Figure 2: Function g{u,(,a) versus u; (a) for a = 0.25 and for ^ = 0.01, ( = 0.1, C = 0.2, ^ = 0.4 
and ^ = 0.6. The value of <l>i = eg/|A| is explicitly indicated for the case = 0.2; (b) the function 
£f(u, is plotted in the region u > j for a = 0.25 and for — 0.1, — 0.2 and = 0.3. The 
values of the edges "I'2 and ^3 are explicitly indicated for the case = 0.2. 



the mean-field equation can be cast in the foUowing form: 



— = g{u,C,a) = U \1 - Ca _ I ' (10) 

where Ca = C/(l + o;)^ and 7=|1 — a|/(l-|-a). As expected, despite the inclusion of an 
additional non-magnetic impurity potential, this result coincides with that obtained in Ref. 1^ 
by diagrammatic expansion in the T-matrix approximation. This equivalence occurs because 
the nature of the electron motion (diffusive or ballistic) has no effect at the homogeneous 
mean-field level. By contrast, soon we will see that the nature of the underlying electron 
dynamics does impact on the nature of the fluctuations around the saddle-point. 

The homogeneous saddle-point solutions are given by the values of u satisfying Eq. ([lO|). 
As in the Abrikosov-Gor'kov theory, a number of features of the DoS can be deduced starting 
from the shape of the function g(u, C, a) when the impurity concentration C and the scattering 
amplitude a are varied. In Figure || the typical dependence of g{u, C, a) in u is shown. The 
parameter 7 defines the second order pole, while for values of C, sufficiently low the function 
shows one extremum in the region u £ [0,7) and the two extrema in the region u € (7, 1], 
that disappear as the value of the impurity concentration increases. As will be clear later, 
the three extrema of g have a clear interpretation in terms of the sharp edges of the DoS. 

As in the Born limit (a = 0, 7 = 1), the maximum of g{u, ^, a) in the interval u £ [0, 7) 
has the meaning of the rescaled energy gap, $1 = eg/|A|, so that the set of the values of 
a and ( for which eg becomes zero defines the gapless region for the superconductor in the 
unitarity limit: 

5'(^,C,a)L=o = C = Co=(l-a)^ 

This means that in the region C > {1 — a)^ the superconductor is gapless (e^ = 0, |A| 0). 
Compared with the Born limit, this result shows that the concentration at which the system 
enters the gapless phase is renormalised by the finite scattering amplitude a. 
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Figure 3: Rescaled energy gap eg/|A| = <l>i, edge of the impurity band $2 and edge of the continuum 
$3 versus (; (a) for a = 0.25, Ci < Co = (1 - af ^ 0.56; (b) for a = 0.6, Ci > Co = 0.16. The dashed 
line is the energy gap in the Born approximation case, that in this case coincides with the edge of 
the continuum. 

The two extrema in the region u G (7, 1], which we denote <1>2 and $3 (see Fig. ^ (b)), 
represent respectively the top of the impurity band and the bottom of the continuum of bulk 
states. The dependence of eg/|A|, <I>2 and $3 on ( is shown in Figure]^ for two representative 
values of a. Defining ^0 and (^i respectively as the value of C at which the energy gap becomes 
zero and the value of ( at which the top of the impurity band touches the bottom of the 
continuum, when Ci < Co (Fig. ^ (a)) the top of the impurity band touches the bottom of the 
continuum before the superconductor becomes gapless, while the opposite situation pertains 
when Co < Ci (Fig. | (b)). 

The assigned definitions of eg/|A|, $2 and $3 become clear when we plot the normalised 
mean-field density of states (i^MF(e))y.s/'^ as a function of the rescaled energy, e/|A|: 



In order to explicitly find (i'MF(e))y,s/'^: we have to solve equation (pO|), expressing u in terms 
of e/|A|. For values of a and C at which all the edges, eg/|A| = $1, $2 and $3 are different 
from zero, equation ( |l0| ) has complex solutions for $1 < e/|A| < $2 and for e/|A| > $3. 
For this reason, in this case, the DoS of a single particle excitation consists of two distinct 
parts: over the energy range $1 < e/| A| < $2 the system exhibits an impurity band while for 
e/l A| > $3 there exists a continuum of bulk states. Within the intervening energy intervals, 
the mean- field DoS is zero. As C increases, the energy gap closes, while the top of the impurity 
band converges on the edge of the continuum. Which happens first depends on the value of 
the magnetic scattering amplitude a. Figure |l| clearly shows the phenomenon of the impurity 
band growth around the energy 

Eb |1-Q;| 
= 7 = -;— , 



lAI 



corresponding to that of the bound state developed around a single isolated magnetic impu- 
rity 
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Figure 4: Integration contours for boson-boson and fermion-fermion fields in the complex 0' = —iO 
plane. 

This concludes our discussion of the solutions of the homogeneous mean-field equations 
and their ramifications on the DoS. However, as mentioned in the introduction, the integrity 
of a quasi-particle energy gap (or gaps) predicted by the mean-field theory does not seem 
tenable. Following the discussion of the Born scattering system jl^, we expect optimal 
fluctuations of the random impurity potential to generate fluctuations in the effective spin 
scattering rate which in turn must lead to the nucleation of localised tail states which soften 
the gap edge. In Ref. it was shown that such localised sub-gap quasi-particle states are 
reflected in spatially inhomogeneous instanton field configurations of the cr-model action. In 
the present case, we expect an analogous situation to persist even in the unitarity scattering 
limit: however, in this case, we must expect that all three of the hard edges $i, $2 and $3 
are softened by the bounce configurations. 

3 Instantons and Tail States 

To explore the structure of the tail state distribution, it is necessary to revisit the saddle- 
point equation (|^) and look for inhomogeneous solutions. For this purpose, it is convenient 
to recast the saddle-point equation (||) in terms of its first integral 



D{Wef - \ A\V{9) = const , 



where 



v{e) 



4U— cos( 



smf 



2a 



In 



2a cos ( 



represents the effective complex potential. The identification of the saddle-point solution 
is simplified by the observation that the homogeneous mean-field configuration (^) satisfies 
the condition 0MF(e < £g) = — vr/2 — icfiuF, with 0mf real (i.e. so that the mean-field DoS, 
(:^MF(e)) v',5 — J^SRcos^MF, vanishes below the energy gap). One can therefore identify a 
'bounce' solution parameterised hy 9 = — 7r/2 — icf), with (j) real, and involving the real 
potential K('/') = V{~tt/2 - i(f>). 

Now integration over the angles 9 = diag(i0BB, 6'ff)bf is constrained to certain con- 
tours p7|: the contour of integration over the boson-boson field 9^^ includes the entire real 
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axis, while the fermion-fermion field iO^F runs along the imaginary axis from to in (see 
Fig. As the saddle-point solution must be accessible by a smooth deformation of the 
integration contour, the bounce solution is accessible only to the boson-boson field, 



i9BB{r) = - 



(r) 



(11) 



i.e. the particular bounce configuration is non-trivial in the superspace. With this parame- 
terisation, the saddle-point equation assumes the form (Vr/^f/))^ -I- Vuicj)) = Vr(0mf), where 



K(</') = 4( — sinh< 



cosh( 



flnll 
2a L I 



- 2acosh20| - l] 



(12) 



Typical shapes of the potential for a fixed values of the scattering amplitude {a = 0.25) 
and the magnetic impurity concentration (C = 0.05) and in two different regions of the energy 
are shown in Fig. ^. A bounce solution with minimum action exists in the regions e < ei — eg 
(Fig. ^ (a)) and £2 < e < £3 (Fig. || (b)), where <I>; = e(/|A|, while outside both the unique 
solution is the homogeneous one, = ^mf- In fact when the energy approaches one of the 
hard edges predicted by the mean-field theory, e ^ ej^, e — )■ or e — )■ £3 , the maximum of 
the potential for (jj = ^mf merges with the minimum of the potential corresponding to the 
lowest value of the action. 

Now, for simplicity, let us first focus on the quasi-one dimensional case. Later we will 
generalise the discussion to the d-dimensional case. The symmetry broken solution ( |ll| ) 
involves the real instanton action 



Sinst = 4.Triy\A\^S^ 



(13) 



where (jjm-v represents the 'classical turning point', Va((/)inv) = V^r(0mf)- As the energy ap- 
proaches one of the three edges (e ^ ej" , e ^ e J or e ^ ej^ ) , the minimum of the potential jl^ ) 
corresponding to the lowest value of the instanton action ( p^ disappears, merging with the 
maximum at the mean-field point (/)mf- Developing the function g{u, C, a) (|l^) around one of 
the extrema ez/|A|, 



4 = (-1)' 



5"«,C,a)| 



1/2 



it is possible to deduce the expansion of the potential Vr(0) in powers of (0 — ^mf) around 
each edge: 



Fh((/.)-K(<^mf) 



(-1)' 



1/2 



(14) 



Here the positive dimensionless coefficients v4;(uq, a) and Bi{uq, a) depend solely on the 
scattering amplitude a and the magnetic impurity concentration ^. 

As the energy approaches one of the hard edges, from below or above depending on 
whether we are dealing with $1 and $3 or $2, the expansion ( p^ permits one to obtain an 
analytic solution for Sa,, 



o 4 A, 

J<t, — 



5/2 



15 Bf 



I I e_f/ 
lAI 



5/4 
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Figure 5: Rescaled potential Vr((j!)) — Vr((/>mf) versus the variable (j)' = (j) — (j)uF for a = 0.25, ( — 0.05; 
(a) e/\A\ = 0.2, e/|A| = 0.3 < eg/\A\ and e/|A| = 0.37622 ~ eg/|A|. The inversion point is indicated 
with (jiinv, while the minimum of the potential is indicated with (pmin, (b) e/|Aj = 0.95 > $2, 
£/|A| = 0.96, e/\A\ = 1, e/|A| = 1.0001 and e/lA| = 1.001 < $3. 



and the bounce solution 



4i{x) - (fit. 



Ai 



Bi cosh x/2rQ 



where the extent of the instanton is set by 



A 



1/2 



-1)' 



lAI 



1/4 



Therefore, on approaching the edge from the 'gapped region' (i.e. e — )■ e^g or e — > e^) the 
size of the 'droplet' ro(e) goes to infinity. Moreover, the coefficient Ai goes to zero at the 
gapless point and analogously A2 and A3 go to zero vifhen the impurity band merges with 
the continuum. 

When in the vicinity of the mean-field gap edges, a generalisation of the quasi one- 
dimensional results to higher dimensions (1 < d < 6) can be developed by dimensional 
analysis. Using the approximation (p^, the bounce configuration is shown to have the 
scaling form 



(r)- 



roie) 



-^/(rAo(e)) 



Substituting this relation in the expression of the instanton action, one finds that: 

d-2 r / , , \ -1 (6-d)/4 



S'inst = 'iadng 



^ 2^(6-d)/2 



(-1) 



lAI 



(15) 



where g = vDL'^ ^ is the bare dimensionless conductance and is a numerical constant^. 
This closes our discussion of the particular saddle-point solution together with the corre- 



5 Specifically, = / a!u[{Vu/)2 + P{^x) - /»(u) 
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Figure 6: Smearing of the gap edges due to optimal configurations of the random impurity potentials. 
The three edges, eg , €2 and €2 become the mobility edges separating the regions of bulk delocalised 
states from the localised tail states. 

spending action. However, to complete the analysis it is necessary to explore the influence 
of fluctuations around the instanton solution. 

3.1 Fluctuations 

Here we only sketch the important aspects of the fluctuation analysis, referring to Ref. 
for a more detailed discussion in the context of the Born scattering limit. Generally, fleld 
fluctuations around the instanton solution can be separated into 'radial' and 'angular' contri- 
butions. The former involve fluctuations of the diagonal elements 9, while the latter describe 
rotations including those Grassmann transformations which mix the bf sector. Both classes 
of fluctuations play an important role. 

Dealing first with the angular fluctuations, supersymmetry breaking of the bounce is 
accompanied by the appearance of a Grassmann zero mode separated by an energy gap 
from higher excitations. This Goldstone mode restores the global supersymmetry of the 
theory. Crucially, this mode ensures that the saddle-point respects the normalisation condi- 
tion (2^[0])vi'^v — 1. Associated with radial fluctuations around the bounce, there exists a 
zero mode due to translational invariance of the solution, and a negative energy mode (c.f. 



Combining these contributions, one obtains the following expression for the local complex 
DoS in the tail regions (e ^ Cg, e ^ £2 and e < £3), 



where Xo(r) represents the eigenfunction for the Grassmann zero mode, y'LS'^/^ is the 
Jacobian associated with the introduction of the collective coordinate and S'inst denotes 
the instanton action (|l5|). Thus, to exponential accuracy, the complex local DoS in the tail 
region becomes non-zero only in the vicinity of the bounce. 

4 Discussion 

This concludes our investigation of the role of quenched classical magnetic impurities on the 
weakly disordered superconducting system. The results above show that the gap structure 
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predicted by the mean-field theory for magnetic impurities in both the unitarity and Born 
scattering limit is untenable. The hard gap edge(s) predicted by the mean-field theory are 
softened by the nuclcation of domains or 'droplets' of localised tail states. The mean-field gap 
edges become mobility edges separating bulk delocalised quasi-particle states from localised 
tail states (see Fig. The latter are induced by optimal fluctuations of the non- magnetic 
random impurity potential which increase the effective spin scattering rate. 

How significant are these results for experiment? The nucleation of sub-gap states will 
lead to the softening of the transition from the metallic to the superconducting phase which 
would be revealed in measurements of the heat capacity close to the bulk Tc. Similarly, since 
the tail state regions are broad on the scale of the coherence length (and broaden close to 
the quasi-particle energy gap), one can expect the existence of sub-gap states to be revealed 
in measurements of the tunneling density of states close to the mean-field gap. However, it 
should be noted that the strength of arising from the sub-gap states will compete with the 
weight arising from dynamical spin fluctuations. The latter will typically give rise to power 
law tails in the sub-gap density of states . 

Finally, our considerations of the density of states focussed largely on energy scales e ~ 
Eg. Here the generating function is dominated by the homogeneous mean- field and bounce 
configurations of the non-linear ci-model action. However, within the gapless phase (i.e. as 
e ^ 0), field fluctuations Q which commute with both cr™ and a™ become massless. These 
fluctuations, which are parameterised by transformations Q = TQmfT~^ where T = IpH ^ 
Isp fS) t and t = 7(^^^)"^7~^, are controlled by a non-linear a-model defined on the manifold 
T e OSp(2|2)/GL(l|l) (symmetry class D in the classification of Ref. The latter 

reflect quantum interference effects in the particle/hole channel and substantially modify the 
low-energy, long-range spectral and transport properties of both the bulk system and the 
low-energy states inside a localised domain. For a comprehensive discussion of their effect, 
we refer to the comprehensive discussions in Refs. |l2|, p^ . 

Acknowledgments: We are grateful to Austen Lamacraft and Robert Moir for valuable 
discussions. 
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